Direct Measurement of intermediate-range Casimir-Polder potentials 
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We present the first direct measurements of Casimir-Polder forces between solid surfaces and 
atomic gases in the transition regime between the electrostatic short-distance and the retarded 
long-distance limit. The experimental method is based on ultracold ground-state Rb atoms that 
are reflected from evanescent wave barriers at the surface of a dielectric glass prism. Our novel 
approach does not require assumptions about the potential shape. The experimental data confirm 
the theoretical prediction in the transition regime. 

PACS numbers: 



Although being the most precisely tested theory in 
physics, quantum electrodynamics (QED) leads only in 
very special cases to measurable forces. One example 
are the well-known Casimir and van der Waals forces 
[3, In addition to their fundamental importance, 
a detaillcd understanding of these forces is crucial for 
testing new fundamental physics at short distances such 
as non-Newtonian gravitational forces 0, HI- Further- 
more, they have important technological implications 
for the development of micromachines with nanoscale 
moving parts [f| @ . Today, Casimir forces between solids 
can be measured with high precision @, H, H|. These 
measurements are all done with objects that are large 
compared to the relevant distances where Casimir forces 
become dominant. Therefore, the underlying theory 
contains the macroscopic properties of the objects, i.e. 
the dielectric functions. Moreover, also the geometry of 
the macroscopic bodies plays an important role due to 
the non-additivity of Casimir forces. A much cleaner 
situation is given when the test object is microscopic. 
This is in good approximation true for a single atom. 
In this case the force by which the atom is attracted 
towards a surface is often referred to as Casimir-Polder 
(CP) force. The shape of the CP potential depends on 
the distance from the surface. While in the limits of 
short and long distances CP-forces can be approximated 
by different power laws [l], [l(| , in the transition regime 
simple analytical expressions do not exist. Here, the 
potential is given by an integral over the frequency- 
dependent dielectric function of the substrate and the 
polarizability of the atom [Til ]. 



In the last two decades many efforts have been 
made to measure CP forces including sophisticated 
approaches such as diffraction and interferom etry of 
cold atom beams at thin transmission gratings jl2l . Il3j 
and quantum reflection of atoms from solid surfaces 
[pfl , [la , [Til \vf\ . In these experiments CP forces were 
studied indirectly by fitting the coefficients of the 
theoretical surface potentials to the measured data. 
Direct methods that have been developed up to now 
can be divided into spectroscopic [H, [T3, H3, [H| and 
kinetic measurements. The latter include reflection of 



FIG. 1: Experimental situation. A repulsive evanescent wave 
potential Few and an attractive surface potential Fsf sum up 
to build a barrier at a distance zb from the surface with height 
E B . 



cold atoms from an atomic mirror [221 ] and surface- 
induced frequency-shifts of oscillating Bose-Einstein 
condensates in magnetic traps (23|. In this article the 
transition regime is probed for the first time by a direct 
model- free measurement. This is done by reflecting 
ultracold atoms from evanescent wave barriers similar 
to previous work [22| . However, here we systematically 
vary the mirror potential and introduce a new data 
analysis which allows for the direct investigation of 
surface potentials at sub-wavelength distances from 
the surface. The new method makes use of a known 
repulsive potential Few( z ) that is added to the unknown 
attractive surface potential Vsf(z) with the goal to 
generate a potential barrier. By varying the strength of 
Few (z) , the height and the position of the barrier can 
be adjusted. The height is measured by reflecting cold 
atoms of a given energy. The position can be determined 
from the derivative of the barrier height with respect 
to the strength of Vew( z )- This last step is the key 
feature for reconstructing the unknown surface potential. 

The experimental situation is shown in Fig. [TJ An 
evanescent wave leaking out from the surface of a trans- 
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parent substrate generates a repulsive dipole potential of 
the form 



Vew = C • P ■ cxp { -2 



(1) 



with a constant Co and the field decay length zq [241 ] . 
The total potential is now the sum of the (attractive) 
unknown surface potential Vsf and the EW Potential 



Hot = Vsf + Vew ■ 



(2) 



If the repulsive potential is strong enough, a potential 
barrier is formed at a distance zb{P) from the surface. 
At the maximum VJ.' ot = 0, which means that 

Vip(z = z B ) = 2C • - • exp (-2—1 . (3) 

ZQ { Z ) 

Furthermore, the height of the barrier is given by 



E B = V SF (z B ) + Co ■ P ■ exp \ -2^ 

zo 



(4) 



In the experiment this barrier height Eb is measured as a 
function of the laser power P. From the data the deriva- 
tive ^jp- is taken and compared to theory as follows. We 
differantiate ^ with respect to the laser power taking 
the inner derivatives into account and substitute eq. ([3]) 
into the result. Note that zb{P) is a function of P. This 
delivers the following equation: 



dE B n { z B 

—— = Co • exp \ -2 — 
dP { zo 



(5) 



The crucial point here is that this derivative depends 
only on the evanescent wave potential and not on the 
surface potential. From this expression the position of 
the barrier can be calculated to be 



(T3\ Z 1 ( 1 dE B 

ZB(P) = -T' ln U'^P 



(0) 



With this knowledge of the barrier height Eb and 
the barrier position zb the unknown surface potential 
Vsf(z.b) can be determined by solving ([4]). 

The experimental task is to measure the height of the 
potential barrier as a function of the laser power. For 
a given laser power this is done by classical reflection 
of an ultracold atom cloud with variable kinetic energy 
from the barrier. The experiment is carried out with a 
setup explained in detail in [251 ] . Here only a short sum- 
mary is given. An ultracold atomic cloud is prepared 
in a Joffc-Pritchard type trap some hundred microme- 
ters below the superpolished surface of a dielectric glass 
prism which is mounted upside down in a vacuum cham- 
ber. Almost pure condensates can be generated with 
some 10 5 atoms. For this experiment however, only very 
cold thermal clouds are prepared in order to avoid ef- 
fects due to the interaction between the atoms at high 



density. The ultracold 87 Rb cloud is held in the mag- 
netic trap at a fixed distance Z\ below the prism surface. 
Then a vertical laser beam (wavelength A = 830 nm) 
that propagates perpendicularly through the prism from 
the top is switched on adiabatically in order to generate 
a dipole trap with radial and axial trapping frequencies 
of Wd,r = 2-7T x 50 Hz and Wd,z < 27r x 0.1 Hz. In this 
combined magnetic / dipole trap a cloud temperature of 
T r~j 100 nK is measured. Now the atoms are accelerated 
towards the surface by a sudden shift of the magnetic 
trapping minimum to a new variable position z% . After a 
waiting time of a quarter of an oscillation period during 
which the atoms accelerate in the shifted trap, the mag- 
netic field is quickly ramped to a constant gradient which 
compensates for the gravitational force. The atoms now 
move towards the surface with a nearly constant velocity 
i>o- It is determined by absorption imaging of the posi- 
tion of the cloud in the first few milliseconds of its motion. 
While the atoms move to the surface they are slightly ac- 
celerated due to residual curvature of the levitation po- 
tential. This effect is taken into account as a correction 
of the measured velocity. The velocity at the surface is 
then given by v 



with u>\ e 



2?r x 4 Hz. 



During the reflection of the atoms the dipole trap guar- 
antees radial confinement. The measured radial FWHM 
width of the atomic cloud is 40 /zm such that the atoms 
are reflected only from the center of the evanescent wave. 
There, the potential barrier reaches its maximum due to 
the Gaussian intensity distribution of the EW laserspot. 
The EW laser is centered around A = 765 nm with a 
spectral width of AA = ±1 nm. The vertical potential 
generated by the dipole trap is very weak and can be 
neglected. After contact with the surface the number 
of reflected atoms is determined by absorption imaging. 
The measurement is repeated for various trap displace- 
ments Az = Z2 — z\ which correspond to different ve- 
locities v. Typical results are shown in the inset of Fig. 
The data points show a gradual decrease from the 
situation where all atoms are reflected (R = 1) to full 
transmission R = 0. The width of this decrease is mainly 
dominated by the Gaussian velocity distribution of the 
atoms corresponding to the temperature of the cloud. 
It is also slightly effected by the imhomogeneous barrier 
height due to the Gaussian transversal intensity profile 
of the EW laser beam. Nonclassical broadening effects 
like quantum reflection and tunneling may play a role for 
the higher laser powers used, where the surface potential 
is steepened by the evanescent wave 26]. However, its 
influence on the center of the decrease and by that on 
the measured barrier height is negligible. Thus the data 
points in Fig. [2] are fit with a model which implements 
only the two classical broadening mechanisms mentioned 
above. The result of the fits for various laser powers is 
plotted in the main part of Fig. [5] (open circles) . For com- 
parison, barrier heights are plotted as derived from the 
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FIG. 2: Inset: Typical measurement of reflectivity R for a 
laser power of P = 134 mW. The fitted barrier height is Eb = 
&b x (3.62 ±0.03) jj,K. Main figure: Dependence of the barrier 
height on the laser power. The open circles are data points 
obtained from curves similar to the ones shown in the inset. 
The solid line close to the data points is a constraint fit to the 
data as explained in the text with the dashed lines limiting the 
95% confidence interval. The three curves (vdW, ret, trans) 
are theoretical expectations for the barrier height assuming a 
nonretarded van der Waals-like potential, a retarded potential 
and the full QED potential. 



theoretical CP potentials valid in the different regimes 

(z > 0, 



(z < I), 



(7) 



with a typical distance I separating both regimes. For 
calculating the potential coefficients Cj the values for 
an ideally conducting surface Cj c are taken from [23] and 
corrected for the dielectric surface. In the case of the 
C4 coefficient the correction is given by a factor <!>(n) 
with refractive index n = 1.512 (2j. In the case of the 
C3 coefficient an often made approximation leads to a 
correction factor of ^tt. Our corrected coefficients are 

C* 3 = 5.8 • 10~ 49 Jm 3 and C* 4 = 5.4 • 10~ 56 Jm 4 . The 
transition length can be estimated from the intersection 
between the retarded and the nonretarded curve to be 
I ~ = 92 nm. In this range we calculate the Casimir- 
Polder potential correctly by solving the full QED for- 
mula (5.39) in 11]. For this purpose the magnetic per- 
meability of glass is /u(w) = ijl = 1 and the atomic polar- 
izability a{oS) follows from [28|]. The dielectric function 
e(u>) is determined from optical data which are available 
for glass in a wide range [291 ] . The barrier height which 
is derived from this correct CP potential fits best to the 
experimental data, although also here a deviation is ob- 
served. This deviation is particularly large for high laser 
powers, where the data exceed the theoretical values. A 
model-free comparison between experiment and theory 
is possible, if the surface potential is extracted from the 
measurement as explained above. The required deriva- 
tive ^jp- w g p[^^)Ip^"^ is taken from a smoothened 



curve that can be obtained by fitting the measured bar- 
rier heights. The fit function must be chosen very care- 
fully in order to stay model-free. A polynomial fit func- 
tion e.g. would implicitly assume a surface potential of 
a certain shape. To maintain generality the values of 
the fit function at the measured laser powers are the pa- 
rameters of the fit. Furthermore, the fit function fulfils 
the following three constraints: (1) the first derivative is 
positive at each point, (2) the second derivative is also 
positive at each point, and (3) the third derivative is 
negative at each point. The physical reason for (1) is 
that an increasing laser power leads to a growing bar- 
rier height. Constraint (2) is equivalent to the fact that 
for increasing laser power the barrier gets closer to the 
surface and constraint (3) means that the rate at which 
the barrier gets closer to the surface decreases. These 
assumptions are valid for any attractive surface poten- 
tial whose attraction is growing with decreasing distance 
from the surface. This is the only assumption on the po- 
tential shape we make. The result of the fit is plotted as 
dots in Fig. [5] [30]. To guide the eye the fit points are 
linearly interpolated. From the fit points the surface po- 
tential is determined by equations Q and ©. The only 
parameters in this calculation are the ones describing the 
evanescent wave: the field decay length z$ is calculated 
from a measured incidence angle of 9 = 43.4° ± 0.1°, and 
a laser wavelength of A = 765 nm to zq — (430 ± 10) nm. 
The proportionality factor Co is calculated from stan- 
dard dipole trap theory .2,4] with a measured beam waist 
of u>o, x — (170 ± 5)/im and wq, y = (227 ± 5)[im. The 
evanescent wave intensity is given (for p-polarized light) 

hv g EW 2 _ 1 . 4 " 2cos W 2 w : t J, „t... 

^ - n (co S (0)2+n2(n<i S in(0)^-1))> Wltil tne rerraC - 

tive index n of the prism. With these parameters the pro- 
portionality constant is C = (1.35 ± 0.05) • 10~ 27 J/W. 



Figure [3] shows the surface potential determined from 
the measured barrier height (dots). Statistical errors are 
due to the spread of the data points in Fig. [2] around 
the fit curve. For the vertical axis they are given by 
a shift of the fit curve to the 95% confidence interval 
boundaries. For the horizontal axis the evaluation of the 
statistical error depends on the error of the gradient of 
the fit curve. An estimation of this error can be given by 
calculating the mean gradient and assuming a maximum 
gradient by a linear interpolation between the lower left 
border of the confidence interval with the upper right 
border. For the minimum gradient the upper left bor- 
der is connected with the lower right border. This gives 
an uncertainty in the gradient S i^jp-)- The uncertainty 
of the barrier position Szb is then given by means of ([6]). 
Systematic errors are due to the uncertainty for zq, Cq, P 
and v. They turn out to be smaller than the systematic 
errors. Taking into account these uncertainties the mea- 
surements agree best with the full QED calculation. The 
retarded and the static potential can be excluded. Other 
experiments show that patch potentials from adsorbed 
Rubidium atoms can play a role in the surface poten- 
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FIG. 3: Measured and theoretical Casimir-Polder potentials: 
in the large figure the theoretical surface potentials are plot- 
ted, i.e. the nonretarded van der Waals potential (vdW) and 
the retarded Casimir-Polder potential (ret). The full theoreti- 
cal curve also valid in the transition regime (trans) approaches 
the retarded curve for large distances and the nonretarded 
curve for small distances. The inset magnifies the colored 
box, in which the measured data points are lying. Statisti- 
cal and systematic errors are indicated by the bright and the 
dark grey shaded area respectively. 

tial [3l| • Such potentials increase the attraction between 
surface and atom. In contrast, our measurements show 
a slightly smaller attraction than theoretically expected. 
Therefore, patch potentials seem to be negligible in this 
setup. This might be due to a permanent exposure of 



the prism surface to the evanescent wave which can ei- 
ther lead to laser induced desorption of atoms [32| or to 
an increased diffusion of atoms on the surface. 

In this article direct measurements of the Casimir- 
Polder force between ground-state Rb atoms and the 
surface of a dielectric glass prism at distances between 
160 nm and 230 nm are presented. A novel method has 
been introduced which is based on a test potential gen- 
erated with an optical evanescent wave at the glass sur- 
face. The measurements do not coincide with the limiting 
formulas valid in the static and in the retarded regime. 
However, they do agree with a correct QED calculation 
and thus confirm the present theory. In addition to the 
mentioned measurement errors the observed small devia- 
tion might be caused by the incomplete knowledge of the 
dielectric function of the used borosilikate glass prism. 
For calculating the theoretical curve the well-known di- 
electric function of SiC>2 glass has been used. However, 
the optical properties of glasses vary depending on the 
exact type of glass [33]. It is therefore possible that the 
theoretical curve slightly deviates from the real situa- 
tion in the experiment. Already a moderate increase in 
the experimental resolution will make it possible to dis- 
cern between such theoretical and experimental errors 
and might reveal new physics that could be hidden in 
the observed deviation. 
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